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ABSTRACT 

The spectral energy distributions for pure-hydrogen (DA) hot white dwarfs can be 
accurately predicted by model atmospheres. This makes it possible to define spec- 
trophotometric calibrators by scaling the theoretical spectral shapes with broad-band 
photometric observations - a strategy successfully exploited for the spectrographs on- 
board the Hubble Space Telescope (HST) using three primary DA standards. Absolute 
fluxes for non-DA secondary standards, introduced to increase the density of calibra- 
tors in the sky, need to be referred to the primary standards, but a far better solution 
would be to employ a network of DA stars scattered throughout the sky. We search for 
blue objects in the sixth data release of the Sloan Digital Sky Survey (SDSS) and fit 
DA model fluxes to identify suitable candidates. Reddening needs to be considered in 
the analysis of the hottest and therefore more distant stars. We propose a list of nine 
pure-hydrogen white dwarfs with absolute fluxes with estimated uncertainties below 3 
%, including four objects with estimated errors < 2 %, as candidates for spectropho- 
tometric standards in the range < g < 18, and provide model-based fluxes scaled to 
match the SDSS broad-band fluxes for each. We apply the same method to the three 
HST DA standards, linking the zero point of their absolute fluxes to ugr magnitudes 
transformed from photometry obtained with the USNO 1-m telescope. For these stars 
we estimate uncertainties of < 1 % in the optical, finding good consistency with the 
fluxes adopted for HST calibration. 



Key v^rords: techniques: spectroscopic 
mental parameters. 
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1 INTRODUCTION 

Good flux standards are hard to get. Ground-based observa- 
tions are limited in accuracy by time-dependent variations 
in transparency of the terrestrial atmosphere. Spaced-based 
fluxes are free from atmospheric distortions, but are more 
diflicult to relate to standard sources in the laboratory. In 
addition, charge-transfer efficiency problems in CCDs be- 
come more acute in space due to bombardment by high- 
energy particles. Despite these issues, the Hubble Space 
Telescope CALSPEC spectrophotometric standards, adver- 
tised to provide absolute fluxes good to about 1-2 % in the 
optical and near-infrared, constitute the preferred reference. 

The CALSPEC absolute fluxes are based on model at- 
mospheres for three hot DA white dwarfs normalized to Lan- 
dolt V-band photometry (Bohlin 1996, 2000; Bohlin, Dick- 
inson & Calzetti 2001). In this calibration, the spectral en- 
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ergy distribution postulated for these stars is computed with 
hydrostatic plane-parallel NLTE pure-hydrogen model at- 
mospheres calculated with the code Tlusty (see 33) • These 
atmospheres consist of pure hydrogen, and are free from con- 
vection. Under the assumption that the physics of DA atmo- 
spheres is well-known, these objects are then promoted to 
space calibration sources (see, e.g., Holberg et al. 1982, 1991; 
Finley, Basri & Bowyer 1990; Bohlin, Colina & Finley 1995; 
Kruk et al. 1999; Sing, Holberg & Dupuis 2002; Liebert, 
Bergeron & Holberg 2005; Dixon et al. 2007). Comparison 
with trigonometric parallaxes available for the brightest DA 
white dwarfs gives empirical support to this procedure (Hol- 
berg, Bergeron & Gianninas 2008). 

Fitting intermediate-resolution spectra of these stars 
provides estimates of the two parameters that deflne their 
atmospheric models, effective temperature and surface grav- 
ity, with a precision better than ~ 1% and 0.1 dex, respec- 
tively. These uncertainties influence the predicted spectral 
energy distribution to a level of < 1 % throughout the opti- 
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cal and near-infrared. The solidity of the theoreticaUy pre- 
dicted relative fluxes can be appraised by comparing calcu- 
lations with independent codes. For an effective temperature 
of 20,000 K, LTE models computed by Koester (e.g., Kilic 
et al. 2007) predict continuum spectra that deviate from 
Tlusty's by only 0.5 % in the optical and near-infrared, and 
the differences reach only 1% at 30,000 K, due in part to de- 
partures from Local Thermodynamical Equilibrium (LTE), 
as illustrates Fig. [T] 

DA white dwarfs are ideal calibration sources for other 
reasons. They are intrinsically faint, and therefore located 
at small distances, which renders their observed fluxes un- 
affected by interstellar extinction. Comparatively, B- and 
A-type stars, also with relatively smooth continua domi- 
nated in the optical and near-infrared by bound-free hydro- 
gen opacity, present several disadvantages. They show weak 
metal lines, are located at larger distances from Earth, and 
some are surrounded by disks (Laureijs et al. 2002), which 
might disturb their spectral shapes in the infrared. Further- 
more, many early type stars rotate fast (see the example of 
Vega; Adelman & Gulliver 1990) and, although this is hard 
to detect in the spectrum when the rotational axis is aligned 
with the line of sight, rotation can distort significantly their 
spectral energy distribution and broad-band colors (Perez 
Hernandez et al. 1999, Aufdenberg et al. 2006). 

F-type subdwarfs such as the SDSS standard (BD -1-17 
4718; Bohlin & Gilhland 2004; Ramirez et al. 2006) are 
abundant in the turn-off of the Galactic halo and thick- 
disk populations, but again, they are located at considerably 
larger distances, and their continua are influenced by both 
hydrogen and H~ photoionisation, with the contribution of 
the latter varying depending on the number of free elec- 
trons available, which makes their spectral shape dependent 
on their exact metal abundance, increasing the complexity 
of the model atmosphere calculations. 

In addition to the three fundamental DA stars that 
set the CALSPEC absolute fluxes, Hubble Space Telescope 
(HST) spectrophotometry links observations of other stars 
(some of which are DA) to the fundamental standards, ex- 
panding the calibration network to about two dozen objects 
throughout the sky. This is important, as three standards 
are not enough for accurate calibrations, in particular for 
ground-based observations where small angular separations 
between target and calibrator are a must. Nonetheless, this 
two-step process may result in a loss of accuracy and, due 
to the miscellaneous nature of the secondary calibrators and 
the resulting large range in distance and extinction, signifi- 
cant variations in quality across the sky. 

The link between the secondary HST fiux standards 
and the three primary ones depends on the calibration of 
the time-dependent response of the HST spectrographs. The 
risk of using secondary standards is nowhere more evident 
than in the recent update of Vega's HST spectrophotome- 
try. After considering charge transfer efficiency corrections 
for the STIS CCD (see Goudfrooij et al. 2006), the slope in 
the Paschen continuum was modified by about 2%. Such a 
change implies a decrease in the inferred effective temper- 
ature from 9550 K to 9400 K (Bohlin 2007), and in fact, 
it solves a long-standing discrepancy between the HST and 
the lUE-INES flux scales (Garci'a-Gil et al. 2005). 

Instead of using secondary standards, based themselves 
of observations that are ultimately tied to the three white 



dwarfs chosen by Bohlin and collaborators, the flux scale 
can be set in a single step using a larger number of well- 
behaved pure-hydrogen hot white dwarfs. But those must 
be first identified. 

While bright fiux standards have been favoured in the 
past, there is a growing need for fainter standards. The 
SDSS camera and spectrographs cannot observe standards 
brighter than about 14 magnitude in their normal mode of 
operation: transfer delay integration CCD operation would 
not be feasible, and very short exposure times would lead to 
large errors in the fiuxes. This situation is likely to become 
more exhacerbated for upcoming instrumentation. 

The Sloan Digital Sky survey is a major source of large 
numbers of DA white dwarfs in the range 14 < V < 21 (Har- 
ris et al. 2003, Kleinman et al. 2004, Eisenstein et al. 2006). 
In this work, we use spectra and photometry in the sixth 
data release (Adelman-McCarthy et al. 2008) to identify hot 
DA stars, deriving a homogeneous set of atmospheric param- 
eters and absolute fluxes for them to propose an expanded 
network of DA primary standards. Section 2 elaborates on 
our analysis procedure. Section 3 describes our sample and 
the results. Section 4 describes the determination of the zero 
points of the absolute flux scale, and in §5 we present our 
standard candidates. Section 6 revisits the determination 
of atmospheric parameters for the HST primary standards, 
and §7 closes the paper with a brief summary. 



2 ANALYSIS 

We focus on DA stars with a surface temperature warmer 
than 20,000 K. This avoids known problems with models for 
DA cooler than about 12,000 K for which unexpectedly large 
surface gravities are derived from the analysis of Balmer 
lines (see, e.g. Bergeron, Saffer & Liebert 1992, Eisenstein et 
al. 2006). Our choice of temperatures avoids the presence of 
convective energy transport in these atmospheres, allowing 
us to safely assume radiative balance. 

We use an unpublished grid of model fluxes for DA 
stars previously calculated by Hubeny. The model atmo- 
spheres contain only hydrogen and are computed, in Non- 
LTE, with the code Tlusty (Hubeny & Lanz 1995, Lanz 
& Hubeny 1995). The code, and the hydrogen model atom 
used are publicly available from the Tlusty web sit^B- We 
considered surface gravities in the range 7 < loggr < 9.5 
(c.g.s. units). Among the most important ingredients, H line 
profiles were calculated with the Stark broadening tables of 
Lemke (1997), and the level dissolution, occupation proba- 
bilities, and the pseudo-continuum opacities were taken from 
Hubeny, Hummer & Lanz (1994). 

Our spectral analysis uses the same code and principles 
outlined in Kilic et al. (2007), but we have added one more 
dimension to the problem in order to account for reddening. 
For each star we determine the trio (Tad, log g,E{B — V)) 
of a model that matches best the observations in the re- 
gion 385 — 540 nm (air wavelengths). The spectra are ei- 
ther normalized by a constant, taken as the median value 
of the flux in the selected spectral window, or by a polyno- 
mial model of the star's pseudocontinuum (see 33] for more 
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Table 1. Parameters derived from SDSS spectra. 
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Figure 1. A comparison between tlie model fluxes for DA wliite dwarfs predicted from the LTE model atmospheres by Koester (private 
communication, black line), and Hubeny (red line). Only the shape of the spectrum is of interest for our purposes, and the fluxes are 
arbitrarily normalized by a constant to match perfectly at log A = 3.67 (4677.352 A). The lower panels show the flux ratios 



details). A wavelength grid equidistant in log A was em- 
ployed. The theoretical fluxes are computed for a discrete 
grid with steps of 2,000 K and 0.5 dex for Tcff and logg, 
respectively, and smoothed to a FWHM (= SX) resolving 
power of X/5\ = 2000, which approximately corresponds to 
that of the SDSS spectra. Reddening is introduced with a 
discrete grid sampling < E{B — V) < 0.095 in steps of 
0.005 mag, using the prescription of Fitzpatrick (1999) with 
R = Av/E{B -V) =3.1. 

The optimization is based on the Nelder-Mead method 
(Nelder & Mead 1965) assisted by Bezier quadratic interpo- 
lation in the grid of synthetic fluxes (see, e.g., Auer 2003). 
Internal uncertainties are estimated, assuming the errors in 
the fluxes are normally distributed, from the diagonal of the 
inverse of the curvature matrix. 



3 SAMPLE SELECTION 

With the goal of identifying solid DA candidates for a net- 
work of flux standards, we scanned the SDSS DR6 photo- 
metric catalog in search for hot objects satisfying {u — g)o < 
0.6 and {g — r)o < —0.2. We also imposed the constraint 
of a moderate estimated extinction towards the sources: 
Ag < 0.827, with Ag from the dust maps of Schlegel, 
Finkbeiner &i Davis (1998). This is approximately equiva- 
lent to E{B — V) < 0.27 mag, but this value would only 



apply to objects outside the Milky Way, or at least beyond 
the extinction layer; we will later restrict our analysis to 
stars with E{B ~ V) < 0.1. To allow a reliable determi- 
nation of the stellar parameters, necessary to assign unique 
theoretical spectral energy distributions, we also limited our 
search to sources in the range 14 < g < 17. Our query to 
the SDSS Catalog Archive Server (CAS; Thakar et al. 2008) 
returned 598 object^. 

Viable DA candidates must show spectra that are well 
matched by DA models. Thus, we downloaded the spectra 
for all the selected targets and analyzed them following the 
prescription in iJH The quality of the fittings was assessed 
from the statistics, and used to select the sources that 
were closely approximated by theoretical fluxes for DA stars. 
Although white dwarfs are very dim, at the faint magnitudes 
covered by the SDSS hot DA are at distances of a few hun- 
dred pc, and extinction becomes significant, at least when 
aiming at modeling flux distributions with high accuracy. 

The radial velocities determined by the SDSS pipeline 
were used to correct the Doppler shifts in the spectra be- 
fore the analysis, and both the model fluxes and the ob- 
served spectra were normalized to have a median flux equal 
to one in the wavelength range used in the analysis. Only 



^ The query returned, in fact, 599 sources, but one file ("spSpec- 
53388-1937-316. flt" ) had an anomalous format. 
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Figure 2. Model fittings to four of tlie SDSS sample stars. The fluxes are in fx units, normalized to have a median of value of one in 
the range 3850—5400 A. Note that the parameters here do not match those in Table[l] as the table gives the average results for two types 
of analyses: those using spectra that preserve continuum information such as those shown here, and those using continuum-corrected 
spectra. 



57 objects which spectra were matched by the model fluxes 
with < 0.63 and with derived atmospheric parameters 
21,000 < Tcff < 85,000K and 7.1 < logg < 9.4, and lim- 
ited extinction, 0.0 ^ E{B — V) < 0.095, were retained for 
further analysis. Fig.[2]shows a few representative examples 
of the fittings in the case where the spectra are normalized 
by their median fluxes in the analysis window. As mentioned 
earlier, we also analyzed the spectra after rectifying the con- 
tinuum shape. 

In Fig. |3] we show the determined atmospheric param- 
eters for the sample. As expected, the derived gravities de- 
crease smoothly for hotter stars (see, e.g. Althaus et al. 
2005). At least in the range 20,000 < T^s < 40,000 K, 
three clear groups are apparent. These likely correspond to 
the three mass clusters found for DA white dwarfs in stud- 
ies using the Palomar Green survey (Liebert et al. 2005), 
and in earlier SDSS data releases (e.g. Kepler et al. 2007, 
DeGennaro et al. 2008). Perhaps due to our strict selec- 
tion criteria for the quality of the fittings, and the resulting 
minute error bars, the three groups appear here more cleanly 
separated than in previous reports. We note that heavy ele- 
ments are more commonly found in the atmospheres of DA 
stars hotter than 50,000 K (Barstow et al. 2003b), making 



such stars statistically more vulnerable to underestimated 
effective temperatures when their hydrogen lines are inter- 
preted with pure-H model atmospheres (Barstow, Hubeny 
& Holberg 1998). The group of low-mass white dwarfs with 
Tcft < 40, 000 K and log g < 7.5 are probably the result of 
stellar evolution in close binary systems, which makes their 
radial velocities likely to vary. 

Kilic et al. (2007) considered both continuum-corrected 
fiuxes and relative fluxes on the determination of atmo- 
spheric parameters for DA white dwarfs from medium reso- 
lution MMT spectra. Their conclusion was that using well- 
calibrated relative fluxes helped considerably to reduce de- 
generacies, and no significant systematic differences in the 
inferred surface gravities were apparent. The continuum- 
corrected spectra are immune to the distortions of interstel- 
lar extinction, and therefore similar surface temperatures 
obtained by the two methods places a limit to the interstel- 
lar absorption that the spectra in our sample have suffered. 

We repeated such an experiment for the particular case 
of the SDSS spectra considered in this work. Reddening 
needs not be considered when analyzing the continuum- 
corrected spectra. However, if reddening is neglected (forcing 
E{B — V) — 0) in the analysis of spectra containing contin- 
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Figure 3. T_eff and logg derived from the analysis of spectra with continuum information such as those illustrated in Fig. [2] The 
uncertainties are internal and do not consider systematic errors in the models or the observations. Note that a non-hnear scale is being 
used for the abscissae for the sake of clarity. 



uum information, a systematic trend is apparent. This is 
illustrated in Fig. UJ where the derived T^s from continuum- 
corrected spectra are typicaUy higher, in particular for the 
hotter stars, in line with the expected signature of inter- 
stellar reddening. When reddening is considered, such trend 
is eflectively removed, and the agreement between the two 
methods greatly improved, as shown in Fig. |4l Systematic 
differences are also present in loggr, for both cases (with and 
without reddening corrections), but these are also reduced 
when reddening is accounted for and have nevertheless a far 
smaller impact on the absolute fluxes. 

The results based on spectra with continuum informa- 
tion are more precise, but are also more exposed to system- 
atic errors in the SDSS flux calibration. Our analysis of such 
spectra, when reddening is considered, returns parameters 
that are on the same scale as those from continuum-rectified 
spectra. Therefore, we adopt as the final parameters (Tcff 
and log g) the weighted average of the two analyses. It is also 
noticeable in Fig. U that the internal random errors derived 
from each method are not always consistent. Conservatively, 
we adopt as 1 — ct uncertainties the sum in quadrature of the 
estimated random uncertainties and the systematic differ- 
ence between the two methods. By fitting the spectra with 
continuum information with the Tcff and logt; fixed to the 



average values and allowing E{B ~V) to vary, we derive an 
estimate for the uncertainty in E(B — V) for each star, which 
is again added in quadrature with the internal uncertainty to 
obtain the final estimate. The final adopted parameters are 
given in Table [T] Typical (median values) uncertainties are 
1.8 %, 0.047 dex, and 0.032 mag in Tch, logg, and E{B-V), 
respectively. These uncertainties in T^s and log g axe slightly 
larger than those inferred by Liebert et al. (2005) for their 
analysis of spectra with typically higher signal-to-noise ra- 
tios but lower resolution. There are three spectra of the star 
SDSS J033133.89+010327.9, and the derived parameters are 
fairly consistent. 

Most of the stars in our sample (exactly 43; see Table 
[l}, are included in the SDSS-DR4 catalog of DA stars of 
Eisenstein et al. (2006). Between 20,000 and 30,000 K, the 
Tcff values of Eisenstein et al. are lower than ours by a few 
percent, and the trend reverses for warmer temperatures, 
with the Eisenstein et al. values being up to 10 % higher 
than ours at about 65,000 K. We find good agreement be- 
tween our surface gravities and those from the SDSS-DR4 
catalog at the lowest temperatures in our sample (~ 20, 000 
K), but growing discrepancies at warmer temperatures, with 
the SDSS-DR4 values being - 0.25 dex higher around 65,000 
K. There are 15 stars included in Table[l]which are also part 
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Figure 4. Comparison between the atmospheric parameters retrieved for our SDSS sample stars using continuum-corrected spectra and 
spectra that preserve the continuum shape. The left-hand panels correspond to the case when reddening is neglected, while the right-hand 
panels correspond to an analysis that accounts for reddening 



in the catalog of DA stars from the Palomar Green survey 
(Liebert at al. 2005). The Liebert et al. surface tempera- 
tures are also a few percent lower than ours in the range 
20, 000 < Tcfi < 30, 000 K, but they are lower by as much as 
10 % at about 60,000 K. Their surface gravities are system- 
atically higher by 0.08 ± 0.02 dex, and it should be noted 
that a similar discrepancy was also found by Liebert et al. 
when they compared their gravities with several other inves- 
tigations (Finley et al. 1997, Marsh et al. 1997, Homeier et 
al. 1998, Koester et al. 2001). In the hot end, our tempera- 
ture scale is intermediate between those of Eisenstein et al. 
and Liebert et al., which difler from each other by ~ 20 %. 



4 ABSOLUTE FLUXES 

Once the atmospheric parameters, surface gravity and effec- 
tive temperature, have been determined for the DA, their 
surface radiative flux is given by the corresponding model 
atmospheres. Only a small part of their spectrum (385 — 540 
nm) is actually used in the parameter determination de- 
scribed above, but the employed grid of model fluxes extends 
to the entire optical region (300 — 700 nm), and therefore we 
determine absolute fluxes for the broader spectral region. 

The quantity of interest is the flux received at the Earth, 
which is related to the stellar surface flux provided by the 
model atmosphere by the angular diameter, or more pre- 



cisely by the square of the ratio of the distance between 
the star and the Earth and the stellar radius. Fairly accu- 
rate estimates for the radii of DA white dwarfs are possible 
based on models, but no direct estimates are available for 
the distances to these faint DA. A direct approach that can 
be followed instead is to use broad-band photometry to scale 
the theoretical flux distributions. Fig. [5] illustrates the span 
of the SDSS ugr bands and the coverage of the theoretical 
spectra. With a typical uncertainty of 0.02 mag in each fllter, 
considering the three together gives an expected uncertainty 
in the absolute scale of the fluxes of about 1%. This figure 
is comparable or smaller than the expected uncertainty in 
the relative fluxes, but it presumes no systematics. 

In practice, we derive synthetic absolute fluxes at the 
stellar surface for the corresponding atmospheric parameters 
and extinction for each DA by interpolation in our grid, 
and determine synthetic photometry in the ugr bands using 
the fllter responses for point sources observed at an airmass 
of 1.3 available from the SDSS web pagefl For clarity, we 
refer to the observed magnitudes with letters in italic {ugr, 
indicating SDSS PSF magnitudes as extracted from the DR6 
CAS) and to those calculated at the stellar surface in bold 
(ugr). For each star, the difference between the magnitudes 
computed for the the photosphere and those measured at 



^ http:/ /www.sdss.org/dr3/instruments/imager/index.html#filters 
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Table 2. SDSS Photometry and zero points. 



Name u g r Su Sg Sr S {g and r) 



SDSS 


J033133 


89+010327 


9 


16 


183 


(0 


022) 


16 


432 


(0 


013) 


16 


876 


(0 


012) 


56 


633 


56 


587 


56 


632 


56 


610 


(0 


023) 


SDSS 


J033133 


89+010327 


9 


16 


183 


(0 


022) 


16 


432 


(0 


013) 


16 


876 


(0 


012) 


56 


717 


56 


651 


56 


675 


56 


663 


(0 


012) 


SDSS 


J033133 


89+010327 


9 


16 


183 


(0 


022) 


16 


432 


(0 


013) 


16 


876 


(0 


012) 


56 


747 


56 


678 


56 


701 


56 


689 


(0 


012) 


SDSS 


J075058 


46+491707 


3 


16 


664 


(0 


015) 


16 


740 


(0 


018) 


17 


117 


(0 


018) 


56 


294 


56 


263 


56 


267 


56 


265 


(0 


002) 


SDSS 


J075106 


52+301726 


4 


15 


649 


(0 


015) 


15 


917 


(0 


020) 


16 


310 


(0 


022) 


56 


351 


56 


265 


56 


169 


56 


217 


(0 


048) 


SDSS 


J081126 


68+053911 


9 


16 


261 


(0 


016) 


16 


416 


(0 


012) 


16 


778 


(0 


010) 


56 


126 


56 


110 


56 


103 


56 


107 


(0 


004) 


SDSS 


J081234 


94+040852 


1 


16 


723 


(0 


016) 


16 


860 


(0 


012) 


17 


302 


(0 


009) 


56 


737 


56 


674 


56 


700 


56 


687 


(0 


013) 


SDSS 


J082346 


15+245345 


7 


15 


273 


(0 


028) 


15 


552 


(0 


025) 


16 


039 


(0 


014) 


55 


984 


55 


911 


55 


916 


55 


914 


(0 


003) 


SDSS 


J084537 


74+065346 


2 


16 


539 


(0 


018) 


16 


637 


(0 


016) 


16 


972 


(0 


014) 


55 


932 


55 


967 


55 


968 


55 


968 


(0 


001) 


SDSS 


J091558 


19+201606 


2 


16 


529 


(0 


013) 


16 


590 


(0 


015) 


17 


026 


(0 


010) 


56 


316 


56 


221 


56 


261 


56 


241 


(0 


020) 


SDSS 


J091601 


87+200758 


1 


16 


674 


(0 


013) 


16 


941 


(0 


015) 


17 


395 


(0 


010) 


57 


160 


57 


130 


57 


123 


57 


126 


(0 


003) 


SDSS 


J092010 


55+045721 


1 


16 


274 


(0 


019) 


16 


673 


(0 


016) 


17 


188 


(0 


010) 


57 


652 


57 


624 


57 


659 


57 


642 


(0 


017) 


SDSS 


J094203 


19+544630 


2 


16 


714 


(0 


013) 


16 


934 


(0 


017) 


17 


398 


(0 


016) 


57 


027 


56 


975 


56 


977 


56 


976 


(0 


001) 


SDSS 


J094940 


37+032425 


5 


16 


454 


(0 


020) 


16 


787 


(0 


016) 


17 


301 


(0 


014) 


57 


656 


57 


571 


57 


595 


57 


583 


(0 


012) 


SDSS 


J095230 


44+114202 


3 


16 


317 


(0 


022) 


16 


472 


(0 


016) 


16 


869 


(0 


019) 


56 


454 


56 


379 


56 


331 


56 


355 


(0 


024) 


SDSS 


J095245 


59+020938 


9 


16 


026 


(0 


015) 


16 


347 


(0 


021) 


16 


832 


(0 


018) 


57 


250 


57 


135 


57 


099 


57 


117 


(0 


018) 


SDSS 


J100222 


50+292755 





15 


677 


(0 


010) 


16 


136 


(0 


Oil) 


16 


669 


(0 


016) 


57 


428 


57 


395 


57 


382 


57 


388 


(0 


007) 


SDSS 


J100543 


92+304744 


7 


15 


967 


(0 


023) 


16 


346 


(0 


018) 


16 


904 


(0 


013) 


57 


651 


57 


544 


57 


558 


57 


551 


(0 


007) 


SDSS 


J100614 


76+441906 


3 


16 


401 


(0 


024) 


16 


824 


(0 


014) 


17 


369 


(0 


014) 


57 


892 


57 


846 


57 


856 


57 


851 


(0 


005) 


SDSS 


J101328 


17+061207 


4 


15 


845 


(0 


018) 


16 


280 


(0 


019) 


16 


754 


(0 


021) 


57 


180 


57 


165 


57 


116 


57 


141 


(0 


024) 


SDSS 


J103743 


48+485720 


8 


15 


275 


(0 


021) 


15 


303 


(0 


021) 


15 


648 


(0 


015) 


54 


475 


54 


501 


54 


533 


54 


517 


(0 


016) 


SDSS 


J 103907 


38+081840 


9 


16 


222 


(0 


007) 


16 


180 


(0 


017) 


16 


595 


(0 


021) 


55 


701 


55 


604 


55 


656 


55 


630 


(0 


026) 


SDSS 


J104332 


62+445329 





16 


271 


(0 


012) 


16 


673 


(0 


014) 


17 


183 


(0 


016) 


57 


644 


57 


591 


57 


573 


57 


582 


(0 


009) 


SDSS 


J104419 


01+405553 





16 


660 


(0 


028) 


16 


832 


(0 


040) 


17 


237 


(0 


017) 


56 


855 


56 


788 


56 


739 


56 


763 


(0 


024) 


SDSS 


J110634 


39+073712 


2 


16 


548 


(0 


Oil) 


16 


878 


(0 


030) 


17 


305 


(0 


025) 


57 


047 


57 


076 


57 


098 


57 


087 


(0 


Oil) 


SDSS 


J114152 


63+253533 


9 


16 


620 


(0 


026) 


16 


990 


(0 


029) 


17 


480 


(0 


016) 


57 


787 


57 


736 


57 


725 


57 


731 


(0 


005) 


SDSS 


J120525 


01+303444 


7 


15 


697 


(0 


024) 


15 


877 


(0 


022) 


16 


300 


(0 


014) 


56 


132 


56 


016 


55 


965 


55 


990 


(0 


026) 


SDSS 


J121205 


11+140801 


8 


16 


810 


(0 


012) 


16 


811 


(0 


020) 


17 


153 


(0 


014) 


56 


186 


56 


139 


56 


132 


56 


136 


(0 


003) 


SDSS 


J121845 


69+264831 


8 


16 


320 


(0 


Oil) 


16 


469 


(0 


012) 


16 


861 


(0 


021) 


56 


272 


56 


237 


56 


211 


56 


224 


(0 


013) 


SDSS 


J122336 


20+412242 


7 


16 


887 


(0 


019) 


16 


926 


(0 


026) 


17 


285 


(0 


020) 


56 


593 


56 


512 


56 


465 


56 


488 


(0 


024) 


SDSS 


J124407 


67+582351 


9 


16 


589 


(0 


017) 


16 


846 


(0 


026) 


17 


300 


(0 


018) 


57 


231 


57 


144 


57 


114 


57 


129 


(0 


015) 


SDSS 


J125217 


15+154443 


2 


14 


339 


(0 


033) 


14 


403 


(0 


029) 


14 


800 


(0 


021) 


54 


294 


54 


189 


54 


150 


54 


169 


(0 


020) 


SDSS 


J130234 


44+101239 





16 


627 


(0 


018) 


16 


988 


(0 


018) 


17 


470 


(0 


024) 


57 


729 


57 


677 


57 


663 


57 


670 


(0 


007) 


SDSS 


J132232 


12+641545 


8 


16 


054 


(0 


017) 


16 


268 


(0 


024) 


16 


658 


(0 


016) 


56 


381 


56 


333 


56 


249 


56 


291 


(0 


042) 


SDSS 


J132434 


39+072525 


3 


16 


417 


(0 


019) 


16 


581 


(0 


018) 


17 


006 


(0 


023) 


56 


586 


56 


513 


56 


491 


56 


502 


(0 


Oil) 


SDSS 


J133207 


33+665453 


4 


16 


669 


(0 


017) 


16 


812 


(0 


026) 


17 


342 


(0 


015) 


56 


869 


56 


767 


56 


844 


56 


806 


(0 


039) 


SDSS 


J133514 


52+505012 


3 


16 


484 


(0 


017) 


16 


781 


(0 


013) 


17 


262 


(0 


013) 


57 


247 


57 


190 


57 


221 


57 


205 


(0 


016) 


SDSS 


J134430 


11+032423 


2 


16 


482 


(0 


015) 


16 


603 


(0 


018) 


17 


005 


(0 


016) 


56 


380 


56 


328 


56 


323 


56 


326 


(0 


002) 


SDSS 


J140327 


76+002119 


6 


16 


410 


(0 


022) 


16 


856 


(0 


017) 


17 


376 


(0 


018) 


57 


981 


57 


962 


57 


963 


57 


963 


(0 


000) 


SDSS 


J142020 


80+521549 


3 


16 


753 


(0 


024) 


16 


855 


(0 


018) 


17 


273 


(0 


015) 


56 


334 


56 


342 


56 


401 


56 


371 


(0 


030) 


SDSS 


J143059 


88+100142 


9 


16 


537 


(0 


022) 


16 


629 


(0 


019) 


17 


044 


(0 


017) 


56 


428 


56 


355 


56 


343 


56 


349 


(0 


006) 


SDSS 


J143105 


74+042215 


6 


16 


808 


(0 


014) 


16 


791 


(0 


020) 


17 


174 


(0 


Oil) 


56 


443 


56 


325 


56 


321 


56 


323 


(0 


002) 


SDSS 


J143315 


92+252853 


1 


16 


810 


(0 


030) 


16 


772 


(0 


021) 


17 


168 


(0 


015) 


56 


509 


56 


373 


56 


351 


56 


362 


(0 


Oil) 


SDSS 


J143443 


25+533521 


2 


15 


921 


(0 


017) 


15 


942 


(0 


021) 


16 


248 


(0 


012) 


55 


390 


55 


354 


55 


294 


55 


324 


(0 


030) 


SDSS 


J144814 


08+282511 


7 


14 


254 


(0 


025) 


14 


398 


(0 


023) 


14 


802 


(0 


014) 


54 


121 


54 


088 


54 


081 


54 


085 


(0 


003) 


SDSS 


J145415 


84+551152 


3 


15 


584 


(0 


013) 


15 


813 


(0 


014) 


16 


275 


(0 


016) 


55 


908 


55 


851 


55 


859 


55 


855 


(0 


004) 


SDSS 


J145600 


81+574150 


8 


15 


951 


(0 


022) 


16 


191 


(0 


020) 


16 


679 


(0 


016) 


56 


569 


56 


475 


56 


476 


56 


476 


(0 


000) 


SDSS 


J150045 


01+621107 


2 


16 


909 


(0 


016) 


16 


961 


(0 


020) 


17 


326 


(0 


020) 


56 


632 


56 


565 


56 


518 


56 


541 


(0 


024) 


SDSS 


J150422 


29+621718 


6 


16 


328 


(0 


017) 


16 


808 


(0 


018) 


17 


277 


(0 


018) 


57 


962 


57 


958 


57 


886 


57 


922 


(0 


036) 


SDSS 


J151421 


26+004752 


8 


15 


483 


(0 


012) 


15 


687 


(0 


012) 


16 


103 


(0 


012) 


55 


603 


55 


579 


55 


583 


55 


581 


(0 


002) 


SDSS 


J152041 


96+495140 


8 


16 


037 


(0 


040) 


16 


268 


(0 


026) 


16 


722 


(0 


034) 


56 


442 


56 


392 


56 


368 


56 


380 


(0 


012) 


SDSS 


,1152839 


42+011300 


1 


16 


168 


(0 


016) 


16 


547 


(0 


019) 


17 


037 


(0 


014) 


57 


195 


57 


200 


57 


266 


57 


233 


(0 


033) 


SDSS 


J160839 


07+074542 


5 


16 


653 


(0 


019) 


16 


690 


(0 


016) 


17 


093 


(0 


016) 


56 


290 


56 


199 


56 


236 


56 


218 


(0 


019) 


SDSS 


J165318 


76+371027 


2 


16 


563 


(0 


019) 


16 


498 


(0 


012) 


16 


843 


(0 


013) 


55 


954 


55 


847 


55 


826 


55 


837 


(0 


Oil) 


SDSS 


J165851 


11+341853 


3 


15 


774 


(0 


Oil) 


16 


173 


(0 


015) 


16 


697 


(0 


012) 


57 


315 


57 


240 


57 


227 


57 


234 


(0 


007) 


SDSS 


J170331 


62+223251 


3 


16 


496 


(0 


016) 


16 


608 


(0 


Oil) 


17 


002 


(0 


013) 


56 


106 


56 


109 


56 


158 


56 


134 


(0 


025) 


SDSS 


J172406 


14+562003 


1 


15 


830 


(0 


021) 


16 


029 


(0 


022) 


16 


419 


(0 


014) 


56 


741 


56 


557 


56 


438 


56 


497 


(0 


059) 


SDSS 


J212412 


14+110415 


7 


16 


835 


(0 


022) 


16 


792 


(0 


013) 


17 


165 


(0 


012) 


56 


424 


56 


295 


56 


276 


56 


286 


(0 


010) 


SDSS 


J214001 


05-075052.2 


15 


756 


(0 


020) 


16 


019 


(0 


012) 


16 


460 


(0 


013) 


56 


212 


56 


174 


56 


164 


56 


169 


(0 


005) 
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2.5x10-' 



1.5x10-' 



1 .0x10-' - 



5.0x10 



3000 4000 5000 6000 7000 

Wavelength (A) 

Figure 5. Model flux for a DA white dwarf with T^f{ = 20, 000 and log g = 7.0. The shapes of the SDSS photometric passbands {ugr) 
are shown with thick lines, arbitrarily scaled. 




the Earth can be derived for each band as Sm = m — m, 
with m(m) replaced by u, g or r (u, g or r). In the absence 
of systematic errors, these three values will be consistent 
and 



Sm 



-2.5 log (5) 



(1) 



where R is the stellar radius, and d the distance to the star. 
(Recall that the effect of interstellar extinction is already 
included in the model magnitudes m). 

Table 2 gives the observed magnitudes and the S val- 
ues for each star. There is an average offset < 5u — Sg >= 
0.0590 ± 0.0056 mag, but the results are more consistent for 
g and r with < 5r - Sg >= -0.0068 ± 0.0051 mag. The ob- 
served u magnitudes should be reduced by 0.04 magnitudes 
to place them on the AB system (Oke & Gunn 1983) , as rec- 
ommended in the SDSS web pages (www.sdss.org; see also 
the discussions in Eisenstein et al. 2006, Ivezic et al. 2007). 
However, after applying this correction, the values for Su we 
derive are still offset from those for g and r by 0.02 mag. 
There is no dramatic variation in the uncertainties in the 
observations for each filter, with median values of 0.0179, 
0.0176, and 0.0150 mag for u, g, and r, respectively, but 
dropping u seems the most sensible option. There are known 
issues with the u filters, such as flatfielding problems, a red 
leak which might affect the observations particularly if there 
are nearby unresolved cooler objects (Adelman-McCarthy 
et al. 2008, Stoughton et al. 2002), and a recently discov- 
ered variation with time of the effective filter transmission 
caused, most likely, by a degradation of the UV coating of 
the «-band CCDs (Abazajian et al. 2009). In addition, even 
if the detected Su — Sg offset is removed, the median la scat- 
ter in the average S is 0.025 mag, but this value is reduced 
to 0.016 mag when only g and r are considered. 



5 THE PROPOSED STANDARDS 

As discussed in the introduction, the uncertainties in the rel- 
ative shape of the spectral fluxes computed for DA stars are 
at the level of < 1 %, but given that we do not know exactly 
the atmospheric parameters and the amount of interstellar 



reddening, additional contributions to the error budget need 
to be considered. 

To evaluate the impact of the uncertainties in the de- 
rived Tcfi, logg and E{B — V) on the optical spectral en- 
ergy distributions assigned to each of the DA, we interpolate 
fluxes after perturbing each parameter by the expected 1 — a 
uncertainty. Because we are normalizing the fluxes of our 
sample stars to match the SDSS g magnitudes, we normal- 
ize the spectra at the weighted average wavelength for the g 
passband, which is 4686 A. Fig. [S] illustrates the results for 
the four stars included in Fig. (2] Depending on the star, the 
uncertainties in Teff or E{B — V) dominate the errors in the 
relative fluxes; logg has always a minor effect. Large errors 
are observed in the cores of the Balmer lines, which are par- 
ticularly sensitive to changes in Toff and surface gravity. We 
combine the three contributions in quadrature to estimate 
the final uncertainties. Although it is possible to determine 
and use the internal covariances among the errors in the pa- 
rameters to estimate the uncertainties in the fluxes, the fact 
is that the errors in the parameters are dominated by the 
contributions from the systematic offsets between the anal- 
yses using continuum-corrected spectra and those in which 
the continuum shape is preserved, and therefore we simply 
assume that the errors are weakly correlated. 

The error budget for each of the DA in our sample com- 
bines the wavelength-dependent uncertainties in the relative 
fluxes, just discussed, with those in the zero point magni- 
tudes -set by the average between the observed and pre- 
dicted fluxes in the g and r bands. Useful flux standards 
should have reasonable uncertainties. We select only DA 
stars with an total uncertainty < 5% at all wavelengths in 
the optical, and < 3% between 4500 < A < 7000 A. There 
are 9 stars in our sample that satisfy these criteria. These 
are compiled in Table [S] where we also provide uncertainty 
estimates at four different wavelengths. We include the ex- 
pected errors in the theoretical spectral energy distributions 
(1 %). The theoretical spectral fluxes at the Earth for these 
stars are given in electronic formal^. 



These are available from jhttp : //hebe ■ as . utexas ■ edu/std/] 
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Table 3. Proposed spectrophotometric standards and expected fractional errors in their fluxes. 



Name 


a 


O"A=3500 


O"A=4500 






SDSS J082346. 15+245345.7 


15.552 


0.016 


0.011 


0.013 


0.018 


SDSS J094203. 19+544630.2 


16.934 


0.015 


0.010 


0.010 


0.010 


SDSS J095230.44+114202.3 


16.472 


0.028 


0.026 


0.026 


0.027 


SDSS J132434.39+072525.3 


16.581 


0.016 


0.015 


0.015 


0.016 


SDSS J143315. 92+252853.1 


16.772 


0.020 


0.015 


0.016 


0.020 


SDSS J145415.84+551152.3 


15.813 


0.014 


0.011 


0.011 


0.011 


SDSS J145600.81+574150.8 


16.191 


0.012 


0.010 


0.010 


0.010 


SDSS J151421. 26+004752.8 


15.687 


0.013 


0.010 


0.012 


0.014 


SDSS J212412. 14+110415.7 


16.792 


0.020 


0.014 


0.016 


0.020 



ATeff 




4000 5000 6000 4000 5000 6000 

Wavelength (A) Wavelength (A) 

Figure 6. Estimated fractional errors in the relative fluxes for 
the four stars in Fig. [2] 

6 THE HST STANDARDS 

As a result of the analysis presented in the previous sections, 
we have identified a number of stars for which the absolute 
fluxes have an expected uncertainty < 2 % throughout the 
optical (see Table 3). An obvious question is how these fare 
in comparison with the HST standards. To answer this ques- 
tion we follow exactly the same steps as for the SDSS stars: 
fitting the STIS spectra normalized to their median in the 
spectral window, or normalized to the continuum, adopt- 
ing averaged Tctt/logg, etc. The original spectra, obtained 
from the CALSPEC web site, have a resolving power about 
1000 in our range of interest (STIS L-mode observations), 
which is half of the resolution of the SDSS spectra consid- 
ered earlier, and this taken into account in the analysis. We 
first assumed that the spectra were at rest, and after find- 
ing the best-fitting models, we determined Doppler shifts 
of +5, -18, and 13 ± 2 km s~^ by cross-correlation, cor- 
rected them, and repeated the optimizations. We refer the 
reader to Bohlin, Dickinson & Calzetti (2001, and references 
therein) for more details about the observations. The derived 
parameters are given in Table 4. Fig. [7] is the equivalent of 
Fig.[2]for the HST standards, and includes the fitting to the 
continuum-rectified spectrum of GD 71. Note that the pa- 
rameters previously used for assigning model fluxes to these 
stars were based on LTE model atmospheres. 

As a reference, we passed the HST fluxes for the three 
standards through the SDSS fllter responses for ugr. The 
results are included in Table 5. Holberg & Bergeron per- 



Table 4. Atmospheric parameters for the three HST standards. 



Name 




logs 


E(B-V) 




(K) 


(dcx) 


(mag) 


G 191 B2B 


61980 ( 514) 


7.555 (0.042) 


0.000 (0.001) 


GD 153 


40401 ( 142) 


7.812 (0.018) 


0.008 (0.001) 


GD 71 


33492 ( 41) 


7.841 (0.017) 


0.004 (0.001) 



formed the same exercise and our results agree with theirs 
to within 0.003 mag, with the exception of the g band for 
GD 71, which they found 0.012 mag brighter. As can be seen 
in the table, the agreement between the derived 5m values 
across the three bands is better than 0.001 mag for G 191 
B2B and GD 71, and only slightly larger (0.005 mag) for 
GD 153, indicating that the good agreement illustrated in 
Fig. [7] extends to the u and r bands. 

To set the zero point of the fluxes we examine first the 
ugr photometry obtained with the Apache Point 0.5-m pho- 
tometric telescope (PT) and transformed to the 2.5-m SDSS 
scale (obtained by S. Kent and D. Tucker; published by Hol- 
berg & Bergeron 2006). The results are included in Table [S] 
We note, as Holberg & Bergeron did, an anomalous g-band 
magnitude for GD 153. After reducing the u magnitudes by 
0.04 mag to place them in the AB system, the zero points 
for G1912B2 and GD 71 agree among the three bands ugr 
to about 0.005 mag. The zero point in GD 153 agrees to 
0.001 mag between the u and r bands (again, after reduc- 
ing u by 0.04 mag), but these disagree with that from the g 
band by 0.05 mag. Comparing the differences between the 
synthetic photometry and the observed magnitudes, the S 
values, with those determined for the HST spectra, we find 
an excellent agreement for G 191 B2B and GD 153 (ignoring 
the g band for the latter). There is, however, a systematic 
offset of about 0.02 mag between the HST absolute fluxes 
for GD 71 and the PT photometry for this star. 

Given these issues, we consider new photometry ob- 
tained using the USNO 1-m telescope. The new observations 
are described in the Appendix, and the mean ugr magni- 
tudes, in the SDSS 2.5-m system, are summarized in Table 
5. We prefer the new photometry and will adopt it in the 
discussion below. The formal uncertainties for the USNO 
photometry {a/N; as given in Table 5) are typically smaller 
than for the PT photometrjjf], due to a larger number of 
measurements available. In addition, while two transforma- 

^ Unlike Holberg & Bergeron (2006), we quote standard errors of 
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Table 5. Photometry and zero points for the three HST standards. The uncertainties quoted correspond 
to the standard error of the mean from the multiple observations available. The PT photometry includes 7 
measurements for G 191-B2B, 1 for GD 153, and 5 in m and 6 in g and r for GD 71. The more recent USNO 
photometry includes 8 measurements in u and 9 in g and r for G 191-B2B, 4 in all bands for GD 153, and 
12 in u and r, and 11 in g for GD 71. 



Name 


u 


9 


r 


5u 


Sg 


Sr 


S {g and r) 


HST spectra 


G 191 B2B 
GD 153 
GD 71 


10.990 
12.661 
12.422 


11.467 
13.065 
12.772 


12.008 
13.578 
13.259 


52.597 
53.711 
53.143 


52.594 
53.704 
53.135 


52.595 
53.713 
53.136 


52.594 (<0.001) 
53.709 (0.005) 
53.136 (<0.001) 






PT photometry (sec Holberg & 


Bergeron 2006) 






G 191 B2B 
GD 153 
GD 71 


11.033 (0.006) 
12.700 

12.438 (0.008) 


11.470 (0.002) 
13.022 

12.752 (0.005) 


12.007 (0.003) 
13.573 

13.241 (0.005) 


52.639 
53.750 
53.159 


52.597 
53.661 
53.115 


52.593 
53.709 
53.117 


52.595 (0.002) 
53.685 (0.024) 
53.116 (0.001) 






USNO 1-m 


photometry (see 


Appendix) 








G 191-B2B 
GD 153 
GD 71 


11.018 (0.004) 
12.700 (0.002) 
12.442 (0.004) 


11.457 (0.002) 
13.064 (0.002) 
12.755 (0.003) 


12.013 (0.003) 
13.592 (0.001) 
13.259 (0.002) 


52.625 
53.750 
53.163 


52.584 
53.704 
53.118 


52.600 
53.727 
53.136 


52.592 (0.008) 
53.715 (0.012) 
53.127 (0.009) 



tions are required to place the PT photometry on the SDSS 
2.5-m ugriz system, only one transformation is necessary for 
the USNO u'g'r'i'z' magnitudes (Tucker et al. 2006). Over- 
all, the scatter between the g and r band scaling factors (5) 
is increased relative to the earlier PT photometry, but more 
in line with expectations at about 1 % (Smith et al. 2002). 
The newer observations bring the p-band magnitude of GD 
153 in line with the scaling factors determined for u and r. 
More importantly, the 2 % offset between the GD 71 HST 
fluxes (i.e. the Johnson V magnitude) and the PT photom- 
etry is now reduced to less then 1%, and therefore within 
the expected variance. 

For the new USNO photometry we find average offsets 
< Su ~ 5g >= 0.0442 ± 0.0017 mag and < 5r - Sg > = 
-1-0.0189 ± 0.0025 mag. Both of these differ by -1-0.02 mag 
from the offsets we determined in §4 for the SDSS DA stars. 
Because those in §4 correspond to measurements made with 
the 2.5-m telescope, we consider them more appropriate for 
the true SDSS system. 

The uncertainties in the relative fluxes of these stars 
due to the uncertainties in their atmospheric parame- 
ters and reddening are smaller than the typical values 
for the SDSS DA stars. This is the result of a much 
higher signal-to-noise ratio of the STIS spectra, and the 
fact that the HST standards are brighter and there- 
fore located closer to the Sun and affected by essen- 
tially no reddening. Nonetheless, four of the stars in Ta- 
ble [5] have similar uncertainties as the HST standards 
(SDSS J151421.26-h004752.8, SDSS J094203.19-f 544630.2 
and SDSS J145600. 81-^574150.8). The uncertainties in the 
absolute fluxes of these stars are dominated by the theoret- 
ical errors, which limit their accuracy at about 1%. 

Our inferred surface gravities for the three DA stan- 
dards are in good agreement with the values reported by 
Barstow et al. (2001, 2003a) from the analysis of Lyman 

the mean in our table, computed from the original table kindly 
provided by D. Tucker. 



and Balmer line profiles, while our effective temperatures 
for GD 71 and GD 153 are a few percent higher than those 
derived in that work. Our temperature for G 191 B2B, how- 
ever, is about 15 % higher than those given by Barstow et 
al. from the analysis of Balmer lines, which range between 
51500 and 54500 kQ. The difference between our analysis 
and the Balmer-based temperatures compiled by Barstow 
et al. can be attributed to the presence of heavy metals in 
the atmosphere of this star (see, e.g. Lanz et al. 1999). Our 
use of a pure-H model leads to an overestimated effective 
temperature (Barstow, Hubeny & Holberg 1998) by about 
6000 K. No metals have been detected in the UV spectrum 
of GD 153 (Barstow et al. 2003b), but we have not found a 
similar study including GD 71. 

In Fig. |5] we confront our absolute fluxes for the these 
stars with those in the CALSPEC library. There is good 
agreement between the spectral shapes, with the differences 
caused mainly by the different effective temperatures used to 
these stars by Bohlin (2000; traceable to Finley, Koester & 
Basri 1997) - our values are higher by 1, 4 and 2 % for G 191- 
B2B, GD 153 ad GD 71, respectively. As explained above, 
our inferred effective temperature for G191-B2B is overes- 
timated by about 6000 K, but it is close to that adopted 
for calibration of the HST fluxes of this star, and hence the 
apparent good agreement shown in Fig. [S] The consistency 
between the HST and our fluxes is best in the red, but over- 
all, the discrepancies are limited to < 1.5 %, in line with our 
expectations. 

Close inspection of the residuals in Fig. [S] ignoring G 
191-B2B, suggests that there is a systematic offset between 
the u (3000-4000 A) and g-band (4000-5500 A) fluxes. The 
HST fluxes in the u band are lower than those presented here 
by about 1 %. The most likely explanation for this feature is 

® We note that the Barstow et al. (2003a) analysis of FUSE spec- 
tra for G 191-B2B leads to higher temperatures (about 57000- 
59000 K), but the differences are believed to be related to diffi- 
culties modeling the Lyman series. 
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.5 



1.0 



0.5 



0.0 



0.5 



0.0 



G 191-B2B 



GD 153 



- Teff = 
logg = 
E(B-V) = 
— \ 



61783 ± 403 
7.572 ± 0.019 
0.000 ± 0.001 
\ \ — 



Teff= 40345 ± 120 
logg= 7.804 ± 0.014 
E(B-V)= 0.008 ± 0.001 
\ \ \ — 



GD 71 





Teff = 
logg = 
E(B-V) = 



33482 ± 50 
7.834 ± 0.010 
0.004 ± 0.001 



Teff= 33503 ± 52 

logg= 7.849 ± 0.010 



4000 4400 4800 

Wavelength (A) 



5200 4000 4400 4800 

Wavelength (A) 



5200 



Figure 7. Model fittings to the three HST primary standards. The fluxes are in fx units, normafized to have a median value of one in 
the range 3850-5400 Ain the top and lower-left panels. The lower-right panel illustrates the model fittings for GD 71 in the case when 
the continuum shape is rectified. Note that the parameters here do not match those in Table [B] as the table gives the average results for 
two types of analyses: those using spectra that preserve continuum information such as those illustrated in the top and the lower-left 
panels, and those using continuum-corrected spectra such as the example for GD 71 in the lower-right panel. The use of pure-H models 
for G 191-B2B is inappropriate as a number of heavy metals have been detected in the UV spectrum of this star. 



the fact that the effective temperatures adopted in the HST 
cahbration for these stars are lower by 2-4 %. 



7 CONCLUSIONS 

In this paper we consider the possibility of building an 
extended network of spectrophotometric standard stars by 
identifying well-behaved pure-hydrogen white dwarfs, spec- 
troscopically determining their atmospheric parameters and 
reddening, and assigning them absolute fluxes based on 
model atmosphere calculations scaled with broad-band pho- 
tometric observations. This procedure was successfully used 
for deflning the three primary flux standards for the HST 
(Finley, Basri & Bowyer 1990; Bohhn, Cohna & Finley 1995; 
Bohlin 2000) , and could be now applied more extensively to 
many other DA stars that are being identified in ongoing 
and planned spectroscopic surveys. 

We put this idea to work using publicly available spectra 
from the Sloan Digital Sky Survey. We identify 59 spectra 
of 57 stars that are nicely reproduced by NLTE DA models, 
and measured their atmospheric parameters. By comparing 



the effective temperatures derived from spectra that have 
been continuum rectified and those which preserve the over- 
aU spectral shape we see the effects of reddening in the hotter 
(and therefore more distant) stars. By including reddening 
in the analysis the two sets of temperatures are brought into 
agreement. 

Our exercise identifies a set of nine new standard can- 
didates with estimated flux uncertainties less than 3% be- 
tween 4500 and 7000 A. We present these stars as poten- 
tially useful standard sources. Their observed spectral en- 
ergy distributions are nicely matched with DA white dwarf 
models across the visible, and although we have not veri- 
fied that their fluxes are stable, all the sources are signif- 
icantly warmer than the ZZ Ceti window (~ 11,000 K), 
and therefore it is unlikely that they show measurable pul- 
sations. Nevertheless, we should exercise caution if trying 
to use them as standards at shorter or, especially, longer 
wavelengths. Some of these stars could have low-mass stel- 
lar or sub-stellar companions which could make their fluxes 
at longer wavelengths to deviate from the predictions from 
the stellar models. There could also be unresolved nearby 
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G 191-B2B 
GD 153 
GD 71 




-0.04^ 
3000 



Figure 8. Ratio between tlie CALSPEC fluxes and those cal- 
culated here using redetermined atmospheric parameters for the 
HST primary standards. The zero point of the flux scale of the 
CALSPEC fluxes is given by Landolt V-band photometry, while 
here it is set by SDSS PT ugr photometry (but see text for an 
exception for GD 153). The two scales are inconsistent by 1.5 % 
for GD 71. 



objects with variable fluxes whose presence went unnoticed 
in the observations analyzed here. 

From the comparison of the scaling factors between syn- 
thetic and observed photometry for the full DA sample we 
find the following offsets: < 5u — 5g >— 0.059 ± 0.006 mag 
and < Sg — Sr >= —0.007 ± 0.005 mag. Assuming the g 
and the r-bands correspond exactly to AB magnitudes (see, 
e.g., the discussion by Eisenstein et al. 2006), and that the 
pure-H models provide accurate descriptions of the spectral 
energy distribution of the DA stars, the derived Su — Sg dif- 
ference could simply correspond to the AB mag correction 
for the u band. However, we note that Holberg & Bergeron 
(2006) found indication of an offset of —0.02 mag between 
the SDSS magnitudes for 107 DA white dwarfs in the g band 
and their fluxes in the Johnson V band. 

We carry out a reanalysis of the STIS spectra for the 
three primary HST standards GD 71, GD 153 and G 191- 
B2B. We use NLTE model atmospheres to determine their 
atmospheric parameters, and find effective temperatures 
that are hotter by 1-4 % compared to those adopted based 
on an LTE analysis. Our predicted fluxes for the three stan- 
dards, scaled using ugriz photometry, agree at a 1 % level 
with the adopted HST fluxes. 

Heavy metals are present in the atmosphere of G 191- 
B2B, and our analysis based on pure-H models results in an 
effective temperature overestimated by about 6000 K. The 
good agreement for GD 153 and GD 71 demonstrates consis- 
tency, but that found for G 191-B2B is simply an artifact and 
highlights that the spectral energy distribution adopted for 
this standard needs to be revised. Our analysis also unveils 
an inconsistency between the photometry obtained for GD 
71 with the Apache Point 0.5-m Photometric Telescope, and 
the zero point of the STIS fluxes for this star, which is based 
on Landolt (1992) V-band photometry, but the discrepancy 
is greatly reduced when considering new u'g'r' photometry 
from the USNO 1-m telescope. Overall, our study supports 
the level of accuracy presumed for the HST standards, but 
we also find room for improvement at the 1 % level. 
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that the g' filter has been determined to have an effective 
wavelength 20 A bluer than that originally quoted). The fil- 
ter transmission data as measured by the Japan Participa- 
tion Group within the SDSS, and the CCD-|-filter response 
curves updated by J. Gunn in 2001 are available on-linqj. We 
refer the reader to Smith et al. 2002 for more details about 
the observations, and to Smith et al. 2007 for a historical 
description of the SDSS photometric system. 

All the individual measurements are given in Table 1X1 
About 9, 4 and 12 observations were obtained for G191 B2B, 
CD 153 and GD 71, respectively. The exposure times were 
(in order u'g'r'i'z') 60, 15, 15, 40, 90 s for G191 B2B, 100, 
25, 25, 60, 130 s for GD 153, and 20, 8, 10, 12, 50 for GD 
71. The average u' g'r'i' z' magnitudes from the USNO 1.0- 
m, included in Table [Xj were converted to the SDSS 2.5- 
m ugriz system using the Equations 2-6 in Tucker et al. 
200qj. We note that the transformations were derived from 
stars cooler than our standards. Rider et al. (2004) found 
that they are still adequate for colors bluer than originally 
defined, although still far from the extremely blue colors of 
the white dwarfs considered in this paper {u — g < —0.3). 
The derived ugr magnitudes are included in Table 5. 



APPENDIX A: USNO U'G'R'I'Z' 
PHOTOMETRY FOR THE HST STANDARDS 

The observations were obtained using the 1.0-m Ritchey- 
Chretien telescope at the U.S. Naval Observatory's Flagstaff 
Station (USNO-FS) during the bright time between October 
2005 and April 2006. The detector was a thinned, UV-AR 
coated, Tektronix TK1024 CCD with a gain of 7.43±0.41 e" 
and a read-noise of 6.0 e~. This CCD is similar to the CCDs 
used in the main SDSS survey camera and the CCD used 
by the 0.5-m photometric monitoring telescope at APO. The 
camera scale of 0.68 arcsec/pixel produced a field of view of 
11.54 arcmin. 

The five filters of the u g'r'i' z' system have effective 

wavelengths of 3540 A, 4750 A, 6222 A, 7632 A, and 9049 7 |http : / /www-star . f nal . gov/ugrlz/Fllter s/response . html] 

A, respectively, at 1.2 airmasses (Fukugita et al. 1996; note * Also given at |http://www.sdss.org/dr7/aIgorithms/jeg-photometrlc-eq-drl.html| 



Faint Spectrophotometric Standards 



Table Al. Photometry for the HST standards in the USNO u'g'r'i'z' bandpasses. 



MJD = JD -240,000,000.5 


u' 


9' 


r' 


i' 


z' 


(days) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


G 191 B2B 


53650 


11.006 


11.517 


12.034 


12.419 


12.739 


53651 


11.024 


11.533 


12.046 


12.421 


12.754 


53654 


11.011 


11.513 


12.019 


12.402 


12.713 


53685 


11.028 


11.525 


12.049 


12.425 


12.739 


53687 


11.035 


11.523 


12.032 


12.422 


12.710 


53702 


11.006 


11.524 


12.035 


12.410 


12.731 


53714 


11.017 


11.511 


12.025 


12.403 


12.735 


53715 


11.017 


11.513 


12.039 


12.419 


12.734 


53740 




11.519 


12.025 


12.406 


12.732 


Mean 


11.018 


11.520 


12.034 


12.414 


12.732 


Std. Dev. 


0.011 


0.007 


0.010 


0.009 


0.013 


Std. Err. 


0.004 


0.002 


0.003 


0.003 


0.005 


GD 153 


53797 


12.695 


13.132 


13.615 


13.974 


14.304 


53832 


12.701 


13.125 


13.611 


13.972 


14.260 


53833 


12.704 


13.121 


13.611 


13.968 


14.300 


53850 


12.700 


13.123 


13.610 


13.982 


14.290 


Mean 


12.700 


13.125 


13.612 


13.974 


14.289 


Std. Dev. 


0.004 


0.005 


0.002 


0.006 


0.020 


Std. Err. 


0.002 


0.002 


0.001 


0.003 


0.010 


GD 71 


53650 


12.422 


12.814 


13.274 


13.646 


13.973 


53650 


12.426 




13.279 




13.987 


53651 


12.451 


12.819 


13.282 


13.656 


13.979 


53654 


12.438 


12.809 


13.271 


13.634 


13.962 
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12.457 


12.819 


13.284 


13.652 


13.965 
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12.423 


12.802 


13.274 


13.644 


13.948 
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12.443 


12.799 
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13.637 


13.967 
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12.820 


13.283 


13.650 


13.966 


53740 


12.444 


12.812 


13.277 


13.639 


13.958 


53742 


12.462 


12.822 


13.295 


13.667 


13.963 


53832 


12.451 


12.828 


13.286 


13.645 


13.976 


53833 


12.445 


12.818 


13.275 


13.651 


13.959 


Mean 


12.442 


12.815 


13.280 


13.647 


13.967 


Std. Dev. 


0.013 


0.009 


0.007 


0.009 


0.011 


Std. Err. 


0.004 


0.003 


0.002 


0.003 


0.003 



